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ABSTRACT 


The increasing use of radar systems as remote sensors of ice 

thickness has revealed a lack of basic data on the microwave properties 

of fresh -water ice. This report describes a program in which the 

complex dielectric constant was measured for a series of ice samples 

taken frcxn the Great Lakes. The measurements were taken at temperatures 
o o o 

of -5 , -10 , and -15 C. It is noted that the ice has considerable 
internal layered structure, and the effects of the layering are examined. 

Values of 3.0 to 3.2 are reported for the real part of the dielectric 
constant, with an error bar of ±0.01. 
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I INTRODUCTION 


The study of naturally occurring ice by airborne radar has been 
continuing since the inception of the early programs in remote sensing 
of earth resources a decade ago. The more recent work by Venier and 
Cross (1972) and Cooper et al. (1974) on lake ice has indicated a need 
for better fundamental data on microwave properties of this type of ice. 
Standard reference works such as Von Hippel (1954) quote values of the 
dielectric constant and loss tangent for ice at a number of isolated 
frequencies and temperatures, but for an understanding of the radar 
signatures of different ice types under naturally occurring conditions, 
more complete data are necessary. The measurements described in this 
report were made on samples of lake ice taken from Lake Erie by NASA 
personnel in the 1973-74 ice season. 
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II SAMPLE DESCRIPTION 


The samples, which were In the form of cylindrical cores 3 inches 
in diameter and up to 36 inches in length, contained several broad cate- 
gories of ice type: 

(1) Clear ice 

(2) Ice with inclusions 0.1 to 1 mm in diameter (milky in appear- 

ance) 

(3) Ice with large Irregular inclusions 0.1 to 0.6 cm in diameter. 

Initially, attempts were made to find long enough sections of each 
sample to be representative of the whole core. This was not possible in 
most cases due to the highly fractured nature of the samples, and an 
alternative method of measuring the properties of the individual layers 
and computing the microwave properties of the whole sample was chosen. 

The NASA samples are described in detail in Table 1. One sample contained 
a small section of ice with verticle columns or pipes of entrapped air, 
but the size of this part of the sample was too small to allow measurement. 

All of the samples contained a substantial number of fractures and most of 
them consisted of a number of distinct layers of different appearance and 
air bubble content. Photographs of sections from four of the samples are 
given in Figures 1 through 4. 

The frequency range employed for this study was 1 (Hlz to 18 GHz. Because 
of the Inherent nature of commercial microwave equipment, it was necessary 
to perform the measurements in three separate ranges — namely, 1 to 8.2 GHz, 

8.2 to 12.4 CUz, and 12.4 to 18 GHz. This in turn made it necessary to 
machine three samples out of each piece of ice to be measured. The low- 
frequency sample was mounted in a 1.625-inch coaxial transmission line, and 
the other samples were machined to fit into an appropriate waveguide section. 


2 



Table I 


DESCRIPTION OF ICE SAMPLES 


NASA 

Sample No. 

Length 

(cm) 

Description 

1 

30 

Clear. Fractured into lengths 
of 6 to 7 cm. 

2 

44 

Inclusion 1 to 8 mm in strata 
10 to 12 cm long. One clear 
section 10 cm long. Fractured. 

3 

47 

Small inclusions, some small 
clear sections. Fractured into 
sections 8 to 12 cm long. 

4 

34 

Mostly clear, occasional verti- 
cal striations. Fractured into 
sections 8 to 12 cm long. 

5 

34.5 

Small inclusions, fractured. 

6 

45 

Mostly clear, fractured into 
lengths of 7 to 8 cm. 

7 

47.5 

Same as Sample No. 6. 

8 

75 

Mostly clear. Bottom section 
crumbled and unusable. 
Fractured . 

9 

42 

Clear, unevenly fractured into 
lengths of 3 to 5 cm. 

10 

78 

Small inclusion (0.1 mm diam- 
eter), 7 cm snow cover. 
Fractured. 7 cm clear sec- 
tion. Occasional strata with 
Inclusions up to 4 mm in diam- 
eter. 
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FIGURE 1 ICE SAMPLE 5 SHOWING SMALL INCLUSIONS 






FSGORE 2 {CE SAMPLE 8 SHOWING CLEAR ICE WITH LAYERS OF LARGE fNeLUSlONS 



fISURE 3 ICE SAMPLE 2 SHOWING HIGH DENSITY OF LARGE INCLUSIONS 
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tA Jf,71 4 

FIGURE 4 ICE EAMPiE 3 SHOWING LAYERS OF INCLUSIONS 
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Ill SAMPLE PREPARATION 


Considerable difficulty was encountered in machining the ice samples 

to fit the various types of transmission line used for the measurements. 

It was found, for example, that the optimum temperature for machining was 
o 

-10 C. Warmer temperatures than this produced local melting as the ice 
was worked, and colder temperatures caused the ice to become excessively 
brittle. Figure 5 shows the coaxial samples being turned down on a lathe 
mounted in an environmental chamber set at -10°C. The dimensional accuracy 
achievable was comparable to that obtained with conventional materials 
(i.e., ±0.001 inch). All samples were made 0.1 m in length. With this 
length of sample, using loss-tangent data from the literature, it was 
estimated that the sample losses would be just above the measurement limit 
of the network analyzer. Longer samples, although preferable, would have 
been difficult to select and fabricate, while shorter samples would not 
give a sufficiently high loss. The coaxial samples were first bored out 
to take the inner conductor of the transmission line, which was then 
frozen into place with the sample symmetrically positioned in the axial 
direction. The outer diameter of the ice was then machined and the outer 
conductor slid into place. Tapered transitions were then added to each 
end to enable the sample holder to be mated with the network analyzer. 

For the frequency region 8.2 to 18 GHz samples were machined to fit 

0.1-m lengths of X- ana P-band waveguide. Initial cutting was done on a 

bandsaw. Further preparation was done with fine sandpaper and a series of 

jigs to ensure a tight fit into the waveguide. It can be seen that for one 

set of measurements over the full frequency range, at least 0.3 m of the 

ice sample had to be of the same type, and had to have no severe fracturing 

present. Following the assembly of the samples in the holders, the samples 

were taken to the analyzer facility and allowed to come to equilibrium at 
o 

-5 C prior to the start of the measuiements. 
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IV MEASUREMENT PROCEDURE 


The measu>'eraent procedure used initially for the Great Lakes ice 
samples is described in detail in the Appendix. In sunnary. the method 
consists of injecting a continuous microwave signal into the sample, and 
measuring the phase and amplitude of both the reflected and transmitted 
waves as the input frequency is swept over the range of interest. The net- 
work analyzer (a computer-controlled H-P 8541A) then calculates the reflec- 
ticm coefficient for an infinite sample, and the loss tangent and complex 
dielectric constant for the material. The procedure could be extended over 
any frequency range between 1 and 18 Qlz in steps of one octave. 

The sample temperature was controlled by two thermoelectric coolers, 
and was continuously displayed during the measurement. Three temperatures 
were chosen; these were -S^"' . -lO"” , and -IS^C. In retrospect, an upper limit 
in temperature closer to the melting point would bave been more appropriate 
in view of typical conditions in the spring on the Great Lakes. 

Some trouble was continually encountered in the precision connectors 
from the sample holder to the analyzer, in that any cooling of these devices 
produced errors in the phase measurements on the sample. This will be dis- 
cussed further under the section on errors. 

The analyzer was limited to approximately 40 frequency points per 
measurement which, in effect, meant that frequency points tended to be 
spaced by 100 to 150 MHz. Fortunately, none of the electrical properties 
of ice appear to be a strong function of frequency in the range covered by 
this investigation, and the spacing was therefore accepted as adequate. To 
start a data run the equipment was calibrated extensively in order to remove 
the propagation effects of the sample holder and its connectors. Further 
checking of the equipment was also achieved by measuring the properties of 
a Teflon sample and comparing the results with known values. 
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After the completion of these formal test procedures the samp e holder 
with the ice in place was connected and allowed to come to equilibrium at 
-5°C. The properties of the ice sample were then measured and the experiment 
repeated for temperatures of -lOf'C and -15*^ C. Because of the time necessary 
to change teaq>erature and come to equilibrium, the measurement of one sample 
would take from 4 to 5 hours, and one day's operation would usurlly result 
in only two samples being completed. Figure 6 shows the sample holder for 
the 1 to 8.2 GHz measurements, and the network analyzer is shown in opera- 
tiem ir. Figure 7. In this figure, the thermoelectric coolers can be seen 
on top of the cold chest. 
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V DATA REDUCTICM 


The majority of the data were reduced using the method described in 
the Appendix. An example of the output from this process appears in Figure 8. 



+.1240 E+05 +1000. 'DIV +.1800 E+05 

FREQUENCY — GHz . ^ _ 

LA-3571 -8 

FIGURE 8 EXAMPLE OF OUTPUT FROM SCATTERING-PARAMETER PROGRAM 
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The options for output format Include plots of loss tangent, losses in dB/n, 
absorption coefficient, permeability, and real and imaginary parts of the 
dielectric constant. Of these, the absorption losses in dB/m and the real 
part of the dielectric constant were chosen as being of the greatest value. 
The permeability of ice was of course unity for all frequencies. In the 
case of data from 1 to 8.2 GHz, an alternative scheme of calculating the 
dielectric properties from the transmitted wave alone was found necessary 
because of the low losses in ice, which allowed multiple reflections and 
standing waves to be set up within the sample. On occasion this allowed 
the ice to act as a quarter-wave plate giving reflection coefficients of 
zero, which in turn caused the conventional program to divide by zero and 
rendered large segments of the data unusable. 

In the alternative program the real part of the dielectric constant 
is given by differentiating the phase of the transmitted signal with respect 
to frequency. This gives the group delay of the wave, and hence the real 
part of e^, as indicated by Eq. (11) in the Appendix. An example of a 
typical phase plot is given in Figure 9 in which it can be seen that the 
slope varies little over the entire frequency range, giving an almost con- 
stant value for c^.. The cyclic variations are caused by the multiple 
reflections within the sample, and have to be removed by smoothing before 
the dielectric constant is calculated. 
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FIGURE 9 EXAMPLE OF PHASE PLOT VERSUS FREQUENCY FOR AN ICE SAMPLE 
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VI RESULTS 


The Ice samples from the Great Lakes contained several different types 
of ice, as explained earlier in this report. All the samples had a layered 
structure, and most included strata of several different ice types. For 
the purpose of this report the ice was classified as clear, milky with small 
inclusions, or clear with large inclusions. Measurements under a microscope 
indicated that the "milky" classification typically contained over 500 
inclusions/cm'^ with diameters of 10 cm, and in most cases also included 
over 9/cm3 with diametersof 5 x 10”^ cm. Using a mixing rule on these num- 
bers, and taking the as 3.2 for pure ice, the bulk dielectric constant 
for the samples would be 3.19 for the small inclusions and 3.12 for the 
samples with both sizes of inclusions. The samples with large inclusions 
typically contained air bubbles with diameters of 0.6 cm at a density of 

O 

0.62/cm**. Using the mixing rule as before, the predicted value for e 

r 

would be 3.05. 

A number of samples of each of these typas were measured, and the 
results are given in Figures 10, 11, and 12. The discontinuities in the 
plots occur at points where the increasing frequency necessitated a new 
sample holder and sample, and correspond to 1 to 8.2 GHz (coaxial trans- 
mission line), 8.2 to 12.4 GHz (X-band waveguide), and 12.4 to 18 GHz 
(P'band waveguide). The error bars on these measurements are larger than 
was hoped, primarily due to the effects of the temperature change on the 
precision cables and connectors and the resonant effects in the low-loss 
ice samples. It is easily seen that in the case of solid ice the real 
part of the dielectric constant is not a strong function of frequency, 
and neither is the dielectric loss, expressed here as attenuation in dB/m. 

An important observation from the standpoint of remote sensing is that the 
ice is virtually los 'ss over the entire range from 1 to 18 GHz, and that 
one should not therefore limit the bandwidth and resolution of airborne 
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FIGURE 10 MICROWAVE PROPERTIES OF LAKE ICE WITH NO INCLUSIONS 


ATTENUATION — dB/m 



FIGURE 11 MICROWAVE PROPERTIES OF LAKE ICE WITH SMALL INCLUSIONS 


ATTENUATION — dB/m 



to 

o 



FREQUENCY — GHz 


LA-3571-12 


FIGURE 12 MICROWAVE PROPERTIES OF LAKE ICE WITH LARGE INCLUSIONS 


ATTENUATION — dB/m 



profiling techniques because of loss considerations In the Ice. The 
losses encountered in the samples that had a significant volume of 
Inclusions are clearly due to a scattering process rather than to di- 
electti.c absorption. 

The magnitude of the losses was for the most part below the noise 

level of the network analyzer, and in all cases was less than 8 dB/m. 

This was a somewhat surprising result and would have no doubt been 

different if we had performed the measurements closer to the melting 

point of the ice. For comparison. Von Hippel (1954) quotes a value of 

o 

1.1 dB/m for pure ice at -12 C, Tables 2 and 3 show the dielectric 
properties of the ten NASA samples constructed from measurements of 
their stratigraphy and the microwave measurements on individual layers 
taken from the samples. One measured sample is not included in these 
results. It was a section from NASA Sample 2, which gave a measured 
value of of 3.65. The measurement was not repeatable du« to deteriora- 
tion of the sample, and the result did not fit the trend for other ice 
measurements and was therefore ignored. 

It was noted that the effect of temperature of the real part of the 

dielectric constant was negligible, and was within the error bars for the 

experiment. This supports the work of Gumming (1952), who maintained that 

temperature was not an important parameter at 10 GHz for ice and snow 

measurements. A trend was noted, however, where e would increase by 
o 

0.002/ C, It is hard to determine whether this was a real trend or whether 
it was due to temperature variations within the connectors and cables. For 
the temperature range under consideration it was in any case a very minor 
effect. 

In the case of dielectric loss, there was a consistent Increase in 
loss as the temperature increased. None of the losses measured were greater 
than 8 dB/m at -5°C. Greater losses would have no doubt been detected at 
temperatures closer to the melting point. 
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Table 2 


VALUES OF Sj. FOR NASA SAMPLES AT C* 


Frequency 

((HIz) 

NASA S 

ample No. 

1 and 9 

2 

3 

4,6, and 7 

5 

8 

10 

1 

3.18 

3.03 

3.06 

3.17 

3.04 

3.13 

3. Cl 

2 

3.18 

3.03 

3.07 

3.17 

3.04 

3.18 

3.04 

3 

3.18 

3.02 

3.07 

3.17 

3.03 

3.14 

3.03 

4 

3.19 

3.02 

3.06 

3.17 

3.03 

3.14 

3.03 

5 

3.18 

3.02 

3.06 

3.17 

3.04 

3.13 

3.04 

6 

3.17 

3.02 

3.06 

3.16 

3.04 

3.12 

3.04 

7 

3.17 

3.02 

3.07 

3.16 

3.05 

3.12 

3.05 

8 

3.17 

3.02 

3.07 

3.16 

3.05 

3.13 

3.05 

9 

3.15 

3.03 

3.06 

3.14 

3.04 

3.11 

3.04 

10 

3.12 

3.03 

3.05 

3.11 

3.03 

3.09 

3.03 

11 

3.14 

3.04 

3.06 

3.13 

3.04 

3.11 

3.04 

12 

3.15 

3.05 

3.08 

3.14 

3.05 

3.12 

3.05 

13 

3.16 

3.05 

3.07 

3.16 

3.06 

3.12 

3.05 

14 

3.18 

3.04 

3.08 

3.17 

3.06 

3.14 

3.06 

13 

3.18 

3.04 

3.08 

3.17 

3.07 

3.14 

3.06 

16 

3.18 

3.04 

3.08 

3.17 

3.07 

3.14 

3.06 

17 

2 18 

3.05 

3.09 

3.17 

3.07 

3.14 

3.06 

18 

3.18 

3.05 

3.09 

3.17 

3.07 

3.14 

3.07 


* 

Values ha«.* been rounded to three significant figures and 
are all subject to a total error of +0.01. For values at 
-IcPc, subtract 0.01 from all readings. For values at 
-lo^C, subtract 0.02. 
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Table 3 

ATTENUATION OF NASA SAMPLES AT -sPc* 
(dB/a) 


Frequency 

(GHz) 

NASA Sample No. 

I and 9 

2 

3 

4,6. and 7 

5 

8 

10 

1 

2.4 

3.6 

2.5 

2.4 

2.1 

2.5 

2.5 

2 

2.4 

3.6 

2.5 

2.4 

2.2 

2.5 

2.6 

3 

2.4 

3.6 

2.5 

2.4 

2.1 

2.5 

2.5 

4 

2.4 

3.5 

2.4 

2.4 

2.0 

2.5 

2.4 

5 

2.4 

3.5 

2.4 

2.4 

2.0 

2.5 

2.4 

6 

2.5 

3.5 

2.4 

2.5 

2.0 

2.5 

2.4 

7 

2.5 

3.6 

2.7 

2.5 

2.4 

2.7 

2.8 

8 

2.6 

3.7 

2.9 

2.6 

2.7 

2.8 

3.0 

9 

2.6 

3.9 

3.4 

2.7 

3.2 

3.0 

3.5 

10 

2.6 

4.3 

3.8 

2.8 

4.0 

3.1 


11 

2.8 

4.8 

3.8 

2.9 

3.8 

3.3 


12 

3.0 

5.3 

3.9 

3.1 

3.7 

3.4 

B 

13 

3.0 

5.3 

3.8 

3.0 

3.5 

3.4 

4.0 

14 

2.9 

5.2 

3.7 

3.0 

3.4 

3.3 

3.9 

IS 

2.9 

5.3 

3.6 

3.0 

3.2 

3.3 

3.7 

16 

3.0 

5.4 

3.5 

3.0 

3.0 

3.3 

3.6 

17 

3.0 

5.6 

3.5 

3.0 

3.0 

3.3 

3.7 

18 

2.9 

5.8 

3.6 

2.9 

3.1 

3.3 

3.8 


These valuer are subject to an error of +1.5 dB/n. 
Values at -ICPc were consistently 1 dB/rc lower than 
the value shown. Attenuation at -IjPc was close to 
the noise level of the Network Analyze and was 
approximately 1.5 dB/m less than the values shown here. 
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VII MEASUREMENT ERRORS 


Two important sources of possible error were present in this 
measurement program — sample preparation, and computational errors. 

In the waveguide measurements (from 8.2 to 18 GHz), the samples were 
accurately fitted and any errors in sample fabrication were due to the 
final facing operation which tended to open up any inclusions near the 
surface and leave nonplanar interfaces. In the case of the coaxial 
measurements (1 to 8.2 GHz) the problem was that the inner conductor 
could produce chipping around the end faces of the sample, and these 
were difficult to "repair" accurately. 

The dominant source of error in the whole program was the fundamental 
problem of measuring phase accurately in the presence of multiple internal 
reflections within the ice sample. Over a bandwidth of 1 GHz, for example, 
an error in the measurement of the phase of the transmitted wave of 1^ pro- 
ducer an error in e^ of 0.003. Under normal conditions, using medium-to- 

high-loss samples, the network analyzer is capable of phase accuracies of 
o 

1 . With ice, however, the low losses allow the generation of standing 

waves within the sample that modify the phase of the emerging wave and 

o 

cause phase measurement errors of i2 . This produces an uncertainty of 
±0.006 in e^. The measurement of losses within the ice sample was affected 
in the same way, but was limited by the noise level of the network analyzer, 
which effectively produced a lower limit of 0.1 dB on all loss measure- 
ments. Since all the ice samples were 0.1 m in length, the lower limit 
was therefore 1 dB/m. An additional uncertainty was introduced by the 
eftects of cooling on the cables and connectors leading to the sample 
area. Experiments on an emply sample holder were conducted to establish 
the magnitude of this effect, and the results lead to a total figure of 
±1.5 dB/m for the errors in the attenuation data. 
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VIII c(m:lusions 


The samples used in this study did not produce any dramatic unexpected 
results. The real part of the dielectric constant varied from approximately 
3.0 to 3.2 and was dependent on the volume of air inclusions present in the 
sample. Temperature did not affect this parameter noticeably. Dieleucric 
loss varied from 2 dB/m to 8 dB/m, and was also primarily a function of 
inclusions rather than temperature. The forward loss was therefore due to 
a scattering mechanism and not to the conductivity of the sample. 

It is reasonable, therefore, to assume that the backscatter in the 
samples would be proportional to the measured attenuation at any given 
frequency for the samples. 

In no case did we find any properties that would explain anomalous 
backscatter as observed by NASA's sidelooking radar, or a sufficiently 
high value for permittivity that could substantiate the value of 4.0 used 
in early flight tests to obtain a closest fit between the ice thickness 
and the short-pulse-radar results. 

It is felt that with further work on new samples, the error bars on 
the microwave data reported here could be reduced substantially by use of 
improved measurement techniques. 
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AutoiiMtic Measurement of Complex Dielectric Constant 
and Permeability at Microwave Frequencies 

WILLIAM B. WEIR, meiibek, ieee 


Ahtrmci Wat Ik* mtrmt «f Ik* c*flw «■< — l*oi*tk tm 



|N’TIODt?CTfON 


T he measurement of complex diHectric cofn^tant 
and complex permeability is required not only for sci- 
entific but also for industrial applications. For ex- 
ample, areas in vhich knowledf;e cd the properties of materials 
at microwave frequencies (a« descnbed by « and m) are applica- 
tions of Qtcrowave heati:'^^, bi<dof(ical effects of microwaxes, 
and nondestructive testing.* 

Numerous measurement iwihods >uitaHe bw different 
ranges of the njutenVai values fd * and m ha\e lieen gi%*en in 
tiw hooks edite«l by Von Hippd (2 1, |3] and in pulffications of 
the American Society for Testing and Materials. It is p<»ssilde. 
howev«^. to rapidly make measurements <»ver the frequence* 
range from 100 MHz to 18 (iHz with a computer-controHed 
network analyzer such as the Hewlett-Packard M»Klel 8540 
series, and by means of appropriate data processing, to de;er- 
mine the complex values of c and n for materials. 

Using the method described in this paper, the complex 
values oi < and $i are determined from measurements made 
diiectly in the frequency domain A somewhat analogous 
method has been de\^oped [7 ] where measurements are made 
in the time domain of the transient response to subnanosecond 
pulses from a dielectric material. With the time-domain mea- 
surement approach, a Fourier transformation is re<)uireu to 
determine c and from the measured transient response. 
Furthermore, with this approach, the frequencies at which t 
and M values are obtained are hand-limited, depending on the 
time response of the pulse and its repetition frequencx • Using 
the -iystem described in thi-< paper, discrete frequencies in 
le>s than 20-kHz step* may l»e selected anywhere within the 
entire 100-MHz to 18-GHz band. 

Actomatk Measi rrufm Svstew 
A computer-controlled network anaUzer is used to mea- 
sure the parameters of a network consisting of a section of 
transmission line containing the sample of material. The trans. 
mission line section may either be waveguide or a TEM trans- 
mission fine. The network is shown schematically in Fig. 1. 
If coaxial-t<f-wavcguide adapters are required, additional 
lengths of waveguide are inserted between them and the 

Manuscript ipcri v cd April 17. 1973, rrviwd July 9. 1973. 

The author is with Stanford Re*rarch Institute, Menlo Park. CaM. 
9402.^ 

> For an examine of a s^-Mem that can he used for r.ondestructive tr*t- 
ift* see |ll. 
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Fig. I . Tunimimon line aection cnataiiBiig dieVctric material- 

sample holder. This is to insure that the higher order evanes- 
cent modes du-? to the coaxial -to- waveguide adapter> are 
significantly attenuated prior to reaching the sample under 
test. 

Under computer contrtd, netw-cu-k analyzer ‘-y>tem calibra- 
tion and measurement are obtain'd at the reference ^anes 
indicated in Fig. I. This is done over a num^wr of predeter- 
mined frequencies at which the complex value> of « and for 
the material are to be determined. The normalized ?-cattering 
parameters (S„'l td the transmission line section containing 
the material are measured at the calibration reference idanes 
at ports .4 and B and corrected for system errors included in 
the calibration data. The measured scattering param^-ters are 
normalized to the characteristic impedance of the trans- 
mission line section. The reflection coefiicienl at the air- 
to-dielectric interface, and the transmission coefficient 
through the material, are found at reference f^anes 1 and 2. 
These coefficients are found directly from the mea^uwl '-can- 
tering parameters after the appropriate phase correction- have 
been applied to account for the shift in the reference planes 
from ports .1 and B to the materia! interfaces (i.e., idanes I 
and 2. Fig. I). 

From the comple.x reffeclion and transmission roefficten:s. 
the computer associated with the network analyzer deter- 
mines the real and imaginarx* parts of the dielectric con-tant 
and perme.ildiity. the loss tangent, and the attenuation |»er 
unit length t.f material. In addition, data taken at •'C' oral fre- 
quencies are to find the average group delay through the 
sample. Average g*x>updelay, in turn, is used to aut<»matiral1y 
resolve phase ambiguities that result when the sample length 
of material is greater than a wavelength in the dielectrir. The 
complex c and m. h*s.s tangent, and attenuation data are auto- 
matically listed by a leleprinier, or rapidly plotted on an 
*V~ Y recorder, or both. 

Data Peocesmsg Techniqves 

The equations u-ed to compute complex f and m from the 
measured irflertion f.5n) and transmission f.S'n) roeffirients 
are presented below together with the equations that reialt « 
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aad n to atteaualioa aad hm Uafeat. As will be seen, an 
infinite number of roots exist in tbe solution for the etfuations 
« and n, but the eonect root is related to tbe lenfth of the 
sample in terms of wavdength within the material. Tbe means 
of determining the correct soiutkm is shown. 

Referring to Fig. I, the propagatioa (actor for a wave 
propagating through the material is defined as [4| 

/> - *-Ti « (I) 

where 

y propagatioa constant, 
a atteauatioa constant, 
d phase constant. 

A time factor of is not explicitly shown in (I). The phase 
constant is equal to 

2w 

d = - (2) 

where X« is tbe transmission tine guide va\‘eiength. Consistent 
with the definition of the propagation factor (H) |. c and §t are 
dehned in terms of their real and imaginar>* parts as follow's: 

« - « («r' — (3) 

^ = (ft/ — (•*) 

The reflection coefficient (T) at the interface between the air- 
hlled transmission line and dielectnc-lilled line when the ma- 
terial samite is infinite in length may be found from the 
measured reflection (5n) and transmission (5n) coeffid'-nts 
for a sample of hnite length (/) 


where 


r = X ± v^x* - 1 

5»i» - ^ I 


X =“ 


25„ 


( 5 ) 

( 6 ) 


The propagation factor P can. in turn, be found from .S,i. 
Ssu and r : 


5ii + .Ssi “ r 

r + 5,,)T 


C) 


The complex dielectric constant and permeability can he de- 
termined from P and F : 


,\* Vx.’ VV L2rl \pjj 

1 + r 

~ /I r 

where 

Xt free space wavelength, 


( 8 ) 

(9) 


‘ In termii of th« air-fi1k*d line characf^istic imprdanee Z* af>d the 
dielectric-fiJIed line diaracte:i#tif impedance Zm. F •(Zd— 

Notr that 5ti • F when / is ir^-itc. 
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Xr cutoff wavelencth of the transmission line section 
« for a TEM line), 

and 



Equation (8) has an infinite number of roots since the 
inucinary part of the logarithm of a complex quantity [F in 
(8)]ts equal to the angle of the complex value plus 2wn, where 
ft is equal to the integer of (f/X^). Equation (8) is ambiguous 
beca u s e the phase of the propagation factor ^ does not change 
when the length of the material ts increased by a multiple of 
wavelength. How'ever. the delay through the material is 
striedy a function of the total length of the material and can 
be ttseu to resolv’e the ambiguity. 

The phase ambiguity is resolved by finding a M^ution for 
c and It from w'hich a value of group delay is computed that 
corresponds to the value determined from measured data at 
two or more frequencies. For this method to work, the discrete 
frequenc>' steps at which measurements are obtained must lie 
small enough so that the phase <.»f the propagation factor (P^ 
changes less than 360* from one measurement frequency to the 
next. With the use of an automatic measurement wy>.ce»n as 
described in ihi-w t>aper. discrete frequency ste|s.. •wmai! cno »gh 
to meet this mtuirement, can easily l»e 'fleeted. The grm.*i 
delay at each frequency may lie computed for each solution 
of c and li a>'^uming that the change^ in t and m are negligi’>h‘ 
over very small increment- of frequency; 



where / is the frequency in hertz and r,.- is the group delav in 
5€cotid* for the *ftb solution of ($) and i9». The mea-ured 
group delay is determined from the slope of the pha-e of the 
propagation factor v-ersus frequency: 


U 


1 

2* ,// 


on 


where # is the phase in radians of F. Accuracy in the de- 
termination of Tg may be increased l»y applying numerical 
differentiation techniques |,>) where the -lope i- computed 
uring data for three or more frequenrie-. The correct rrsit, 
n»k, is found when 


^ 0. 

Once the correct value- of ' and fi have U'en found at each 
of the measurement frequencies, the las- tangent and attenua- 
tion per unit length may l»e determined. In genera). h»r lM>th 
an electrically and magnetically lo— y material, a h»— tangent 
is defined by* 

tan 4 * — (12) 


where 


a/ 




Mr *r 


tf 


* Ncrtr for the when 1 *^"“®- tan 4»#r" *»'. 
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W'm^ t. 4Mtesrtrtc co«H«a»R« assd prr*«c*fetlttif ^ m tltt^ 

i«| flelatiw cneeitam CN IUHbtiir«' per- 

memiil^y ^1. 



' ««4j-a«£s mt* 


*■ 

* ^ n tt s? 

- S>2 





Pit- I- ^irlpctrjr aod s*erm«abife»:y <sl E<x‘«io«-b SP- 

trw p^meal«lrty 


m 4 

a/^ « i»f 4- 

TW atlensiatMin m p^tr le'tigt.N is found frum th« 

t«'»n 

«»— Ivi + tan**- 1}«>«. (tJ| 

h 

MR4»tmf,» Resi’LTS 

Ttif rt^uli-. i4 tr.i-ajMfing I hr * and # aj" Trfl»« and 

En>rr«.>« and Cutninj; K«T<*^»rl> SF>*.S aw <!«»»■« ia F»g« 2 
thr«ut;li 4 .'t5«'-itrrmrn(» «f l*»ef* (nat«rtai> wrw obtained 
with a 'anijde in'Mi>rt«l in .V-!*and i Wi-S2) waa'ejsitide and at 
an amJ<jVi»t trm}»rfata«‘ >«t Z2*C. The Tertmi ‘#3»|de lengtli 
was I in, the Erfisstirh samjdr Irngih in, 

Tefttn ha» a nf),»tii.el>' low t.’ slAt |&|l in the wtefoware 
regifn and has a very h*»- lt»ss *{ X band -fso *<0.OOtM {jjl. 
ft i' a nuftmaRnrlir malrfial and hrn« Its reJative jierojt: 



»»«»««* «»* 
fmnufMry — isw* 

w 

¥ne. a. Mearrs^ attrmiatim and !«» tan«m5 ts# Errm^srh SP^I.y ss 
tlie X t* ■ ,aj regioa, la} Attemfcstnm itt}.. ianai tanin*m itan Int |, 

aldli'v <|t,i U I The measared da»a shown in Ft|> i eorre'jtund 
ete«ty to these known pro|iertie>i erf T*R«n. 

Some notireaWr e-ariaiK»ns is the 9. I* to 9 *.<,;}!* fre- 
tjoeney laate are indifaittl in Fig. i, htmever. froHj the aver- 
age values <>f t and it Furtherraore, value* imisstng in Fig, 2} 
obtaioetl at three (renMeocies wererrroneoU'. The*e erfi'Mfs are 
attributed t« the fact ‘has the nelworfe analyter 'ignai «t,iurce 
«a* not fr«i|t»ency vtaWiiaeif .»t the time measttrement- were 
obtained «»n the TeB«n material,* VVithoij' a fmjuesev- 
stahilijted -igoa! «>urre, the fre<.|u«ti€*es at the time o( tsea- 
Vtiretnent d« mit eeaclly rejteat those at the time of ralihra- 
ti«i, which ytrtidace* emit* in the (»ha*e of the nte,»'Ur*?d scat- 
tering («araaH*ters. The atxuracy of the t aSues of coraplex » 
and m depends -igmlitaniK on the accuracy of the rot «i«(iresl 
scattering paramcict phase data. 

The ficc<»«»rl> isF-.*,* «* an ahsorfier tyi:* matcftal which 
resenatts. at 5„s t',}lr. The fta'ir comjeMtkro of the material i* 
siik-««e rutifier {tresuniafily It-aded with l>»>iy magneiw- parti* 
cIks When iMcfced with a metallic surface, the material i- re- 
portetl to t« effective is al>o>ri»tng a normally incident wavr 
at S.S ('.Hr. although some aljoirption m-<;ur' at other (re 
tjuenciirs and other angles- of incidence as well. The re<|uire- 
roent <*f a ntefallic surface backing sugge-st* that the absorlirr 
is a magm-sicitlly tos-y material which ahsoflo the high rur- 
wais esisting along the wetatlic surface e.vcifed f»y the inri- 
dent wave. That the material is laaguetically hts-y is con- 
firmed by Fig. 3 which indicates a value of p," greater than 
ssroover X band. Fig, ia!>o indicates that the materia! is not 
electrically !ct««-y fr/'Svff ttver the Itandi but do*-* h.«sf ,i ret.t 
tively high dielectric coBslaut: !«,’2el,f7 owr tht l.jnt)! 
Furthermore, p,' is greater than 1,9. varying from a'cmi 1 •> 
to 1,4 <*ver the Ijund ,‘\,ttenuati«n a«»d loss t.tngetit’* .««■ -h->wn 
in Fig. 4 for the Eccreiorb material. ,''ftiefi»,»N..rt i.uic- fr- m 
about 24 la 40 dUs'em and h,*ss tangent from .jls.ut ti * to ij 9 
aver X hand. 

Ec«««ioiffe SF'S.I 

» t<»l» t^ifSfem <«am if 
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Complex « and |i data have been presented (or two com^ 
mercialiy a\atiaUe dieiectric materials as examples o( results 
obtained using the automatic measurement system described 
in this paper. This system, however, has also been (ound use- 
ful in determining the dielectric properties o( rock, s<mI. and 
other aonsundard materials. 
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